A quantum description of bubble growth in a superheated fluid 
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We discuss a quantum description of bubble growth in a superheated liquid Helium by addressing 
the problem of operator ordering ambiguities that arise due to the presence of position dependent 
mass (PDM) in this system. Using a supersymmetric quantum mechanics formalism along with the 
Weyl quantization rule, we are able to identify specific operator orderings for this problem. This is 
a general method which should be applicable to other PDM systems. 



PACS numbers: 



o 

(N 

C 
3 



43 ■ 
i 

cr 



> 

(N 

q 

o 

o> 
o 



I. INTRODUCTION 

The study of bubble nucleation and growth is a broad 
field of study with rich physics [![. Some of the more re- 
cent examples include experimental and theoretical stud- 
ies of superfluid Helium in a superheated state 0, and the 
theoretical study of nucleation of magnetic bubbles Q. In 
this paper we look at the problem of the quantum de- 
scription of bubble growth in superheated liquid Helium. 
This is just one example of many systems whose classical 
nature is well known but whose quantum representation 
is not. The fact that not all functions can be quantized 
uniquely into operators in the Hilbert space has been 
known for some time. The Groenewald-van Hove the- 
orem has explicitly shown that Dirac's program of as- 
signing operators to all functions by having the Poisson 
brackets change into commutators works only for poly- 
nomials of quadratic form or lessQ. 

In the case of superheated liquid Helium, the difficulty 
in quantization appears in the form of a position de- 
pendent mass (PDM), whereby the kinetic term in the 
classical Hamiltonian is space dependent. In fact, many 
systems in nature have the PDM structure; for instance, 
free carriers, such as electrons, in semiconductors of non- 
uniform chemical composition are often described by a 
PDM Hamiltonian 0, @. Other examples of PDM ap- 
pear in various instances of nuclear many-body prob- 
lems, quantum dots, and nano-mechanical systems just 
to name a few0. As a result there have been a number of 
attempts to address the issue of operator ordering ambi- 
guities due to PDM, using the Gallileian transformation 
and other methods [8]. 

We aim in the present paper to provide insight into 
several existing techniques in quantization, their conver- 
gence and possible future usage. In particular, we show 
in this paper how to apply the Supersymmetric quantum 
mechanics (SUSY QM) formalism!! which have been 
previously applied to solve simpler PDM problems p^l- 
Il2j . The main difference is that, unlike the usual exam- 
ples of SUSY QM, our Hamiltonian contains a specific 
potential which is not easily factorizable; the method pre- 
sented in this paper shows how one may apply SUSY QM 
to a more general situation. 

On the other hand, we apply to our system a more for- 



mal quantization procedure which takes the Heisenberg 
commutation relations as primary, the Weyl quantization 
otherwise referred to as the "Weyl Transform" 0, [HI, [13] ■ 
This allows one to quantize a PDM Hamiltonian but does 
not provide explicit operator ordering. By combining the 
two approaches namely the SUSY QM method and the 
Weyl Transform, we are able to identify specifically a cou- 
ple of correct operator orderings out of infinitely many 
possible orderings. 

The paper is divided as follows: In Section 2, we give 
a brief description of nucleation in superheated liquid 
Helium, which makes the problem of PDM obvious. In 
Section 3, we describe the details of using the SUSY QM 
formalism, to address the PDM problem. In Section 4, 
we introduce Weyl quantization and use it to obtain an 
operator form for our classical Hamiltonian. This result- 
ing Weyl quantum Hamiltonian is converted to a point 
mass representation to compare with the Hamiltonian ob- 
tained through SUSY methods. In Section 5 we study the 
effective potential using realistic parameters, and provide 
a conclusion. 



II. BUBBLE NUCLEATION IN SUPERHEATED 
LIQUID HELIUM 

A superheated (or supercooled) fluid[l5| is typically 
metastable since it cannot make a direct, uniform tran- 
sition to the stable phase throughout its volume. The 
transition to the stable phase occurs via the nucleation 
of a droplet of stable fluid, such that once the nucleus 
reaches a critical size it grows quickly, ultimately con- 
verting the entire liquid from metastable to stable phase. 
Nucleation can occur due to fluctuations in pressure and 
volume within a liquid, whereby small quantities of a new 
phase ("bubbles") form within a homogeneous phase. If 
the liquid were stable these bubbles would disappear; but 
superheated liquid Helium is metastable, so these bub- 
bles can grow and eventually transform the system from 
the old liquid phase onto the new vapor phase. 

This growth process is driven by pressure gradients 
across the surface of the bubble. A recent study [f6j] pro- 
vided a comprehensive classical description of such a bub- 
ble nucleation process that takes into account the fluctua- 
tions in both the radius and the pressure. The Helmholtz 
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free energy corresponding to this process is: 

AF = 2 -lR^ l{Pl -P v f- { P-P v )^- 
+A-kcjR 2 



(1) 



where R is the radius of the bubble and P the pressure 
outside of the bubble. P,; = P + 2a /R is the pressure 
inside the bubble where a is the surface tension, P v is 
the internal equilibrium pressure of the bubble, and Kj 
is the isothermal compressibility of the stable phase at 
the equilibrium pressure. The first term does not affect 
the rate of nucleation since it describes the pressure fluc- 
tuations which average out to a constant. The second 
and third terms of Eq. (|TJ) provide the barrier, with re- 
spect to the radius, over which the nucleus must pass 
in order to expand and fill the volume with the stable 
phase. The critical radius for vapor nucleation is then 
R c = 2a/(P v - P). 

To establish the (classical) Hamiltonian for bubble nu- 
cleation, we note that the kinetic energy of a growing 
nucleus is given by[l[ 



E k 



where the variable mass is given by: 



!/(/,') = 4;t( ■'■ — ~T ) PlR 3 



(2) 



(3) 



with pl being the density of the liquid and p v <C pl 
is the vapor density. The potential term based on Eq. 
([T]) leads to the Hamiltonian for the bubble nucleation 
problem [l7j 



Hcias 



P 



2M Q x 3 



+ U Q x 2 (l - x), 



where x = R/R c is the radius of the bubble scaled to the 
critical radius and Uq = 4iraR 2 and Mq = AttplR^,- Re- 
turning to our particular case in question, we note that 
if we were to quantize our classical PDM Hamiltonian 
naively by replacing x and p by the corresponding posi- 
tion and momentum operators x and p, an infinite set of 
operator forms are possible. For example, ^ — > p^kp, 



h \ ■ M> 2 - h -+ 7/'T/'t to list just 

some of the simplest cases; we show how to identify the 
correct specific operator ordering below. 



III. APPLICATION OF SUSY TO PDM OF 
NON-SYMMETRIC POTENTIAL 

While our potential is not super-symmetric, the super- 
symmetric quantum mechanics (SUSY QM) formalism 
permits the re-writing of the Hamiltonian in terms of 
the familiar creation-annihilation operators by introduc- 
ing an effective potential ("superpotential"). The kinetic 



part of of a general PDM Hamiltonian (with h = 1, and 
without any confining potential) can be written 

H = -lm{x) a ^-m{x) 2b ^m{x) a , (5) 
2 ax ax 

where a (and b) can take any value as long as a + b = — s . 
One can then write the creation-annihilation operators 
A± as 



1 . sh d 



A- = — m (x) b —m(x) a + W a (x) 



(6) 



~m{x) a ^m(x) h + W a (x) (7) 



where W a {x) is known as superpotential and the corre- 
sponding harmonic Hamiltonian is given by 



Hf = At A* = T± + Vt 



(8) 



where 



= -\m{x) a ^-m(x) 2b ■^ z m{x) a (9) 



dx 



r<-> = --m(x) b ^-m(xf a ^-rn(x) b (fO) 
2 dx dx 

are the kinetic terms and are the corresponding po- 
tential terms that effectively play the role of the quadratic 
confining potential of a standard harmonic oscillator. 
The form of W a is fixed by demanding that A^ obey 
the Heisenberg algebra i.e. [A~ , A£] — 1. This condition 
leads to 



f 



W a = ^ / V2Mx)dx+ } 



4a + 1 / f 



\j2m{x) 



,(H) 



(4) which then implies 



4a +1 1 



4 y/m(x) \ ^Jm(x) 



(4a +1) 2 / d 1 



dx yjm(x) 



(12) 



We emphasize here that the procedure is independent of 
the form of the mass - the mass m(x) is completely gen- 
eral. The essence of the SUSY procedure is therefore the 
introduction of an extra potential V a which now contains 
all the effect of operator ordering via the parameter a. 
For our specific PDM corresponding to the case of super- 
heated liquid Helium, m(x) — Mqx 3 so that: 



W a 



3(4o+l) f 



V32M x 5 5 



^V2Mi 



x J 



and 



V, 



21 + 48a - f 44a 2 2 , r - 1 



32Af x 5 



H MnX* =p ■ - 

25 2 



(13) 



(14) 
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It is noted that when a = b = — i 



W a = lf V^&jdx, and U± = ^l^\ (15) 

which are identical, respectively, to the equivalent super- 
potential and potential obtained in the classical formal- 
ism using Poisson brackets in place of commutators [loj . 

To solve our problem, we need to write the Hamil- 
tonian corresponding to our situation in terms of the 
operators A^. Writing out the "momentum" part of 
A^ using our PDM m(x) — Mox 3 , one can show that 
At + A~ = ^y / 2M x 5 which does not depend on the 
parameter a i.e independent of ordering. Using this re- 
sult, one is able to write down the potentials and 
UqX 2 (1 — x) in terms of A+ + A~ in a straightforward 
manner. In terms of these operators, the Hamiltonian 
for our system is given by subtracting V Q + from A+ A~ 
and building the Uqx 2 (1 — x) potential by remembering 
that x is unitless in the original Hamiltonian - this is 
equivalent in this case to multiplying our i by a factor 
(8M0/25) 1 / 5 to get: 



Hr 



-\{At + A-f + U«{A + a +A-f* 



l-(At+A-) 2 / 5 



(16) 



The Hamiltonian of Eq. (fTo]) is our PDM Hamiltonian 
when A^ is defined as in Eqs. (J5J17J) with the superpoten- 
tial W a of Eq. (|13l) . Effectively, we have re- written our 
Hamiltonian solely in terms of the new raising and low- 
ering operators A+ and A~. We note again that in Eq. 
(IT6| . the kinetic energy term is still implicitly with an 
indeterminate operator ordering i.e. A+ A~ — = T 1 -^ 
ofEq. ©. 

We note that the operators A+ and A~ function as 
raising and lowering operators and are constructed such 
that they abide to the Heisenberg algebra, [A~ , A+\ = 1, 
in the same way as the well-known raising and lower- 
ing operators for the standard quantum harmonic os- 
cillator (QHO). Guided by the definitions of the stan- 
dard QHO, we therefore create an ansatz whereby we 
define new "position" and "momentum" variables: x' = 

J^(A+ + A~) and p> = iJ^{A+ - A~). x> and 

p' then play the role of position and momentum oper- 
ators with the usual commutation relation, [x',p'] = i. 
It is noted that the commutation relation for x' and p' 
implies p' = i-4- via the Stone- von Neumann theorem [18| 

i, there is a representa- 
. a d 



[a,b] 



which states that given 

tion such that a\i/j) = x\ip) and p\ip) = In these 

definitions, uj and Mq correspond to the mode energy and 
the mass of the bubble. Scaling our Hamiltonian by mul- 
tiplying throughout by hw and changing the variable to 
z = x' / R c for a later comparison (note that the variable 
z in this section relates to x' . Here z is not a coordinate 



transform of the original position variable x as done in 
the next section), we are led to the following Hamilto- 
nian by a straightforward substitution of x 1 and p' into 
Eq. P]>: 



H eff. - - 



2M0M dz 2 



+ V a {z) + V sys (z) 



where 



and 



Va{z) 



2MnR 2 



21 + . 



144a 



21 



V S y S {z) = U Z 4/5 



1 



100z 2 



7 2/5 



(17) 



(18) 



(19) 



Here Uq is modified to take into account the scaling of 
variables. The effective potential, V a {z)+V sys (z), is fairly 
complicated compared to the original potential; in partic- 
ular it contains a singularity at z = — such is the "price" 
of removing the PDM from the kinetic term. The energy 
spectrum of the effective Hamiltonian of Eq. (|17l) is iden- 
tical to that of the original PDM Hamiltonian (up to a 
constant factor). On the other hand, because the posi- 
tion variable z is now in a different space, the eigenstate 
of Eq. (fT7|) is not the eigenstate of the original Hamilto- 
nian, although there are ways to map the eigenstate of 
Eq. (fTT|) back onto the eigenstate of the original PDM 
Hamiltonian. Once the parameter a is determined, the 
unique operator ordering is decided; we now suggest a 
way to select the value of this parameter. 

Wc note here that the steps involved to derive Eq. p7|) 
is that of "picking a particular basis" for the Hamiltonian 
of Eq. (|16p rather than a coordinate transformation from 
one position coordinate system to another position coor- 
dinate system. In particular, by choosing ansatz based on 
the well-known quantum harmonic oscillator, the result- 
ing Hamiltonian is assured of being of unit measure and 
Hcrmitian. These requirements of unit measure and Hcr- 
miticity is explained in more detail in the next section, 
where a coordinate transformation is directly involved. 



IV. COMPARISON WITH WEYL 
QUANTIZATION 

There exist in the literature several ways to formally 
convert classical function into operator form. Weyl quan- 
tization otherwise referred to as the "Weyl Transform" is 
a well-known way of doing so. There are many qualities 
of the Weyl transform that make it an optimum way to 
perform this conversion [J, Il3l Il4j . The Weyl transform 
is defined as: 



W[FT</>](x) = 



1 

2^ 



[FT](a,P)e lhaP,2 e l0x 



X(j)(x + ha)Q(a, (3)dad/3. 



(20) 
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Here 



[FT](a,/3) 



J T{p, x) exp[ 



-i(pa + x/3)]dpdx (21) 



denotes the Fourier transform of the classical function, 
T{x,p) and 4>(x) is the wave function upon which the 
Weyl operator acts. The weighting factor, il(a,j3), char- 
acterizes the type of quantization to be utilized. Weyl 
quantization, in particular, is defined when fi(a, /3) = 1. 
Therefore, the Weyl transform of the kinetic part of our 
classical Hamiltonian is 



W 



F 



I 2M Q x 3 



p 



Using the result 

r,2 ' 



ha)dad/3. 



F 



P 



2M x 3 



F 



P 



2n5"{a)] 



-sgn(3 



(22) 



(23) 



the Weyl transform for our case is 



F 



P 



(~-2\ 



2M m 



Id 2 3d 3 

x a ox ox x° 



(24) 



We note here that the Hermiticity of the Hamiltonian 
must be ensured. One can show that, given a general 
Hamiltonian of the form H = A{x)-j^ + B(x)-§^ + C(x) 
the Hermiticity condition, 



ip*H(j>(x)dx = / Hip(x) 4>{x)dx (25) 



imposes the following conditions after integrating vari- 
ous terms by parts: (i) A*{x) = A{x) i.e. A(x) is real; 
(ii) A'(x) = Re[B(x)}; (iii) Im[C(a;)] = ±lm[B'(x)] where 
Re[-] and Im[-] denote real and imaginary parts respec- 
tively. This general result shows that the Weyl transform 
such as the result of Eq. (|2"4")l preserves the Hermiticity 
of the Hamiltonian. 

In order to compare this result with the Hamiltonian 
obtained in SUSY, we begin with Eq. ([24| and note 
that the associated kinetic term in the corresponding La- 
grangian is 



L K (x,x,t) = ^M Q R 2 c x 3 x 2 . 



(26) 



We introduce a coordinate transformation x — > cz a 
where c and a are real numbers transforming Lk{x, x, t) 
to L K (z,z,t) = \M R 2 c c'a 2 z 5a - 2 z 2 . In order for the 
l/x 3 factor in the first term to disappear, one needs to 



choose a = | and c = (§) 2 ^ 5 i.e 



2/5 



so that 



the mass term becomes a constant, m(x) — ¥ 1 and Eq. 
(l26l) is transformed to 



L K (z,z,t) 



1 



-M R\ 



2 z 2 . 



(27) 



With this substitution, Eq. (|24|) is transformed into 



H K (z) 



2M R 2 



dz 



3_d_ 

5z dz 



12 



25z 2 



(28) 



Once the coordinate transform is carried out, the mea- 
sure (or the Jacobian of integration | | ) associated with 
the Hamiltonian of Eq. (128)1 is no longer unity i.e. coordi- 
nate transformation does not preserve the inner product 
of the Hilbert space, which is problematic. In general, 
unit measure can be restored through the comparison of 

the inner products between two Hamiltonians H and H 

where H is the Hamiltonian for the inner-product space 
with unit measure and H corresponds to that with non- 
unit measure /j,(x) so as to preserve the Hilbert space: 



i>*(z)H4>(z)fi(z)dz= 4>*(z)H<t>(z)dz 



(29) 



Noting that f(z) — f (z) / yj n(z) where f(z) = <j)(z) or 
ip(z), Eq. ([2U]) implies the relationship between the two 



Hamiltonians as (H* 



general unit measure Hamiltonian H as 



d 2 



d 



H = A(z)-^+B(z)—+C(z), 
oz z oz 



Writing a 



(30) 



and the non-unit measure Hamiltonian as one without 
the tildes, the above relation between H and H implies 



A(z) 
B(z) 



A{z), 
B{z)~ 



A(z) 



H'(z) 



C(z) = A(z) 
-B(z) 



3^'(z) 2 ln"(z) 
2 



4 n(z) 2 
1 f/(z) 



2 n{z) 



fj,(z) 
C(z). 



(31) 
(32) 



(33) 



We note that these results are general results for any 
Hamiltonian of the form Eq. (f3T)]) ; it is not a property of 
a specific Hamiltonian, or any transforms involving dif- 
ferential operators, but is rather a self-consistent result 
following from preserving the inner product, and ensur- 
ing Hermiticity. A corollary of these results is that for 
a non unit measure Hamiltonian, the Hermiticity condi- 
tion implies B(z) — A'(z) + ^J^ A{z) i.e. the Hamilto- 
nian of Eq. (|28p is still Hermitian. This relation also 
means that a point-mass, unit measure Hamiltonian will 
always have B(z) — due to Hermiticity condition i.e. 
the first order derivative term always vanishes, which is 
what we normally see in "ordinary" Hamiltonians. In 



our case, since x 



2/5 



the measure is calculated 

to be n(z) = || = {■§z) i ^ ■ The Hamiltonian with unit 
measure emanating from Eq. (|28p is then 



H(z) 



2M a R 2 



dz 2 



100z 2 



+ t/ z 4/5 (l - z 2 '*) 



2/5N 



(34) 
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Equating this result of Eq. (|34| with the Hamiltonian 
obtained through SUSY QM, Eq. (fTT). we require 



(a) 



(b) 



9 _ 21 + 48a - 144a 2 
100z 2 ~ 100z 2 



(35) 



which implies that specific operator orderings that cor- 
responds to the Weyl quantization rule are given by 
a — > — g or a— Therefore, the quantum Hamil- 
tonian for the bubble nucleation problem that obeys the 
Weyl quantization rule is either 



H = - 



or 



H = - 



1 d 1 d 1 



2Afni? 2 



a; 1 / 2 da; a; 2 da; a; 1 / 2 



da; a; 6 da; 



+ U x 2 (l-x) (36) 



*7 z 2 (l 



(37) 



We have therefore narrowed down the possible operator 
orderings to two (out of infinity). This is the central 
result of our paper. 

We note here that it seems on first sight that one could 
avoid the operator ordering ambiguity from the outset by 
transforming the classical Hamiltonian to a point mass 
form first, and then quantizing it. However, the proper 
procedure is quantization and then coordinate transform 
as we have done in this paper. This is because one loses 
more information than necessary otherwise - the extra 
information that a Weyl Transform gives us as a result of 
having PDM is lost by carrying out the coordinate trans- 
formation first. More specifically, one would not obtain 
the term proportional to 1/z 2 if the point mass transform 
was done first. This would also be inconsistent with the 
SUSY procedure of Section 3 which independently gave 
the term proportional to 1/z 2 . 

Finally, we note the deeper connection between the 
methods of Section 3 and this section. The fact that the 
relatively simple procedure based on SUSY QM formal- 
ism matches the result of Section 4 that involves Weyl 
transform, coordinate transform for finding a point mass 
representation and re-scaling to unit measure is, indeed, 
a demonstration of mathematical self-consistency. At a 
deeper level, the match is not surprising since the re- 
lationship between x and x' of Section 3 is indeed the 
same as the coordinate transform of Section 4, namely 
x' = (2/5)a; 5 / 2 . However, the method of Section 3 is 
somewhat more efficient than that of Section 4 since it 
bypasses the need to find the point mass representation 
and to restore unit measure. This efficiency is possi- 
ble since the quantization takes place not in the position 
space but in the space of raising and lowering operators; 
picking a useful position representation related to well 
known quantum harmonic oscillator then ensures Her- 
miticity and unit measure owing to the Stone-von Neu- 
mann theorem. 




(c) 



> 
+ 

[D 

> 




FIG. 1: (a) Vsys(z) for pressure P = 0.8P„ (Dashed line) and 
P — 0.95P V (Solid line), (b) Same as (a) but for V a (z) instead, 
(c) Va(z)+V 3 y S (z) for the same values of pressure, (d) V a (z) + 
Vsys(z) for several values of pressure P £ [0.8P V , 0.95P V ]. In 
this figure, the position z is in scaled units (in units of R c ) 
and the potential is in eV. 



V. CONCLUDING REMARKS 

The goal of this paper was in finding the correct op- 
erator ordering for the bubble nucleation problem. A 
full quantum mechanical treatment of the bubble nucle- 
ation process that addresses tunneling through the bar- 
rier will be presented elsewhere [lijj. In this section, be- 
fore we conclude the paper, we attempt to gain some 
physical insight into this problem by briefly examining 
the effective Hamiltonian of Eq. ([T7]) [and Eq. ([31|] with 
realistic parameters. We consider physically reasonable 
parameters for this system so that we can get an order 
of magnitude estimate for the parameters. The typical 
values for superfluid Helium at temperature T = AK 
is cr = 0.12 x lCr 3 iV/77i, P v = 8.1445 x l0 4 N/m 2 , 
Pl = lA0kg/m 3 . For zero applied pressure, the classical 
critical radius is R c = 29.5 x lCT 10 m. At T = AK, the 
thermal de Broglie wavelength of a single Helium atom 



and hence ther- 



is A = h/\/2nmkT sa 4.36 x 10" 
mal momentum of pxh = h/A = 1.52 x lO^^kgms^ 1 ^S> 
^/UqMq over all range of applied pressure P. Typically 
a single bubble contains an order of 100 Helium atoms, 
or momentum of roughly lOpxh- 

With these realistic numbers, we have calculated 
and plotted the shape of V a (z) 

V S ys(z) — Uoz 4/>5 (l — z 2 / 5 ) in Fig. [T] for various values 
of applied pressure P/P v . The non-confining shape of 
the potential makes it clear why the numerical solution 
to the eigenvalue problem is difficult and non-convergent. 
The changing shape of V a (z) near z = due to changes in 
pressure and subsequently R c is also shown in Fig. [T] On 
the other hand, the potential for larger z corresponding 



K 9 a „A 
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to V ays . (z) looks simpler. It is clear from Fig. Q] demon- 
strating unbounded potential that the energy spectrum 
is expected to be continuous, and the bubble growth is 
similar to a ball rolling downhill with the "gradient" de- 
pendent on the applied pressure. We note that had we 
not transformed the Hamiltonian and left the PDM in 
the kinetic energy term we would not be able to come up 
with such an intuitive understanding of the system. 

To conclude, we have studied the bubble nucleation 
problem in a superheated liquid Helium and provided a 
quantum mechanical description. Although the system 



seems rather complicated at the outset and furthermore 
involves an exotic Hamiltonian that diverges at the ori- 
gin due to the specific form of the PDM, we were able to 
extract useful information through the methods of SUSY 
QM and the Weyl transform. We have presented a con- 
vergence between two methods by converting to a point 
mass description of the quantum Hamiltonian. Further 
analysis will need to be performed in regards to com- 
parison to the microscopic description of the problem at 
hand, and these will be presented elsewhere jl9j. 
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